Contact interaction of two bodies can be modeled using the penalty function approach while its accuracy and robustness are directly associated with the geometry of contact bodies. Particularly, in the research fields of ...